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Abstract 

For the one-dimensional Hubbard model with Aharonov-Bohm-type magnetic flux, we 
study the relation between its symmetry and the number of Bethe states. First we show 
the existence of solutions for Lieb-Wu equations with an arbitrary number of up-spins 
and one down-spin, and exactly count the number of the Bethe states. The results are 
consistent with Takahashi's string hypothesis if the system has the so(4) symmetry. With 
the Aharonov-Bohm-type magnetic flux, however, the number of Bethe states increases 
and the standard string hypothesis does not hold. In fact, the so(4) symmetry reduces 
to the direct sum of charge-u(l) and spin-s^(2) symmetry through the change of AB- 
flux strength. Next, extending Kirillov's approach [12,13], we derive two combinatorial 
formulas from the relation among the characters of so(4)- or (u(l) © sZ(2))-modules. 
One formula reproduces Essler-Korepin-Schoutens' combinatorial formula for counting 
the number of Bethe states in the so(4)-case. From the exact analysis of the Lieb-Wu 
equations, we find that another formula corresponds to the spin-sZ(2) case. 

1 Introduction 



Low-dimensional physics has attracted a great interest of theoretical and experimental physi- 
cists for almost a half century [16, 22, 23]. Among them, the Bethe ansatz method, which was 
originally developed as a non-perturbative method for diagonalizing one-dimensional spin-^ 
isotropic Heisenberg spin chain [1], opened a new realm of mathematical physics. Roughly 
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speaking, Bethe's work consists of three parts: i) Bethe states are introduced and Bethe equa- 
tions, which are sufficient conditions for the Bethe states being eigenstates, are derived; ii) the 
existence of solutions for the Bethe equations are discussed in simple cases, and the general 
forms of solutions are conjectured, which are called string hypothesis; iii) under the string 
hypothesis, the formula for counting the number of the Bethe states is provided, which leads 
to the combinatorial completeness of Bethe ansatz. At present the Bethe ansatz method 
is applied to various kinds of one-dimensional spin chains and strongly correlated electron 
systems [16,22,30]. 

In Bethe's work, the combinatorial formula for counting the number of Bethe states pos- 
sesses a wealth of mathematical implications. In the case of the isotropic Heisenberg spin 
chain, the Bethe states constructed by finite-valued solutions of Bethe equations do not pro- 
duce all the eigenstates. In fact the system has sl(2) symmetry and the Bethe states are 
sZ(2)-highest [6]. The eigenstates other than highest weight vectors, i.e., s/(2)-descendant 
states, are constructed by applying the lowering operator to the Bethe states. Mathemati- 
cally, the number of Bethe states is interpreted as the multiplicity of irreducible components 
in the tensor products of two-dimensional highest weight s/(2)-modules. Bethe's formula is 
also extended to the generalized Heisenberg spin chains with higher spins or sl(n) symmetry, 
for which the powerful tools such as Q-systems are introduced [12, 13]. 

The application of Bethe ansatz method to the one-dimensional Hubbard model was given 
by Lieb and Wu [21]. The Bethe equations for the Hubbard model are often called Lieb- 
Wu equations. Takahashi's string hypothesis asserts that, in the thermodynamic limit, the 
solutions of Lieb-Wu equations are approximated by string solutions, and the number of Bethe 
states is estimated under the hypothesis [28, 29]. The Hubbard model with even sites has so(4) 
symmetry [7,32,33]. Essler, Korepin and Schoutens proved that, when the system has the 
so(4) symmetry, the Bethe states are so(4)-highest [5]. They also showed in a combinatorial 
way that all the eigenstates are obtained by taking the so(4)-descendants of the Bethe states 
into account [3]. On the other hand the completeness of Bethe ansatz for the system with 
odd sites, which has just sl(2) symmetry related to the spin degrees of freedom, has not been 
discussed. 

In this article, we study the Bethe states in the one-dimensional Hubbard model. In 
particular, we deal with the system on a ring with Aharonov-Bohm-type magnetic flux. The 
system has so(4) symmetry only at special values of the AB-flux strength and the so(4) 
symmetry reduces to spin-sZ(2) symmetry for other values. More precisely, for a generic 
value of the AB-flux strength, the so(4) symmetry breaks into the direct sum of the charge- 
u(l) and the spin-sZ(2) symmetry. Varying the AB-flux strength, we investigate solutions 
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of Lieb-Wu equations. Here we recall that all the enumeration of Bethe states we have 
mentioned above are based on the string hypothesis. However, the violation of the string 
hypothesis is numerically observed: for the spin-i isotropic Heisenberg spin chain, some of 
the string solutions reduce to real solutions when the number of sites is large, which is called 
redistribution phenomenon [4, 8, 9]. Thus, without making any approximation, we discuss the 
existence of Bethe ansatz solutions with the AB-flux. In particular, we show the existence 
of solutions of the Lieb-Wu equations with an arbitrary number of up-spins and one down- 
spin. Here we employ a graphical approach [2, 24]. We exactly count the number of solutions 
in the case and verify that the enumeration with the string hypothesis is correct only in 
the so(4)-case. We find that the Lieb-Wu equations for the system with only the spin-s£(2) 
symmetry have more solutions than those in the so(4)-case. Next we study the combinatorial 
completeness of Bethe ansatz. We obtain the relation among the characters of so(4)-modules 
through the power series identities similar to Kirillov's [12, 13], which gives a new proof of 
Essler-Korepin-Schoutens' combinatorial completeness of Bethe ansatz [3]. We also introduce 
a new combinatorial formula derived from the relation among the characters of (it(l)©s!(2))- 
modules. The formula suggests the combinatorial completeness of Bethe ansatz for the system 
only with the spin-s£(2) symmetry, which has not been discussed in the literature. 

The Bethe ansatz solutions with the AB-flux should be quite important in the low- 
dimensional physics of the one-dimensional Hubbard model. As pointed out by Kohn, the 
low frequency conductivity is directly related to the shift of the energy levels due to twisted 
boundary conditions [15]. Sharpening Kohn's argument on electron systems in any dimen- 
sions, Shastry and Sutherland discussed effective mass of the one-dimensional Hubbard model 
through the twisted boundary conditions [25] . Furthermore, Kawakami and Yang obtained an 
explicit expression for the effective mass of the electric conductivity for the one-dimensional 
Hubbard model [10, 11]. For the Bethe ansatz solutions with the twisted boundary conditions, 
there are other aspects such as persistent current associated with the AB-flux. 

The article is organized as follows: in Section 2, we review the symmetry of the one- 
dimensional Hubbard model and the Bethe ansatz method. We also describe how to construct 
eigenstates other than the Bethe states. In Section 3 we prove the existence of solutions for 
the Lieb-Wu equations with one down-spin and exactly count the number of Bethe states 
in varying the strength of AB-flux. In Section 4 we study the combinatorial formulas for 
counting the Bethe states in terms of the characters of so(4)- and (u(l) © s/(2))-modules. 
The final section is devoted to summary and concluding remarks. 



3 



2 Bethe ansatz method 

2.1 Hubbard model and Lieb-Wu equations 

We introduce the Hubbard model on an L-site ring with Aharonov-Bohm-type magnetic flux 
<I>. Let c| and q s , (i G Z/LZ,s G {T?l}) be the creation and annihilation operators of 
electrons satisfying {cj S ,Cjt} = {cJ S ' c j*^ = ^ an< ^ i c * s ' c j*^ = < % < ^ s *' anc ^ define the number 
operators by nj s := c| s q s . We consider the Fock space V of electrons with the vacuum 
state |0) {dimV = 4 L ). The one-dimensional Hubbard model is described by the following 
Hamiltonian acting on V: 



where we assume </> := 5>/L G R/27rlR and [/ > 0. It is clear that = H^^n- Furthermore 
has the same energy spectra as those of H^^n. Hence we often restrict the region of <j) 
to ^ (f) < ^jj in what follows. 

Let N be the number of electrons and M that of down-spins. We assume ^ 2M ^ N ^ 
L. Let {ki\i = l, 2, . . . , N}, (Re(fcj) G M/27rlR) denote a set of wavenumbers of iV electrons and 
{A Q |a = 1,2, ... , M} that of rapidities of M down-spins. Given a set of spin configuration 
{si\i = 1,2,..., iV} with N — M up-spins and M down-spins, the Bethe state with {ki, X a } 
has the following form: 

IM;s)t,M= E ^M( s ; a )4i >S i42,«2-"4 JV) » JV |0>- ( 2 - 2 ) 

{l<Xi<L} 

The coefficients i/Jk,x{x;s) in (2.2) are explicitly given in [31]. The Bethe states (2.2) are 
eigenstates of the Hamiltonian (2.1) if {ki,X a } satisfy the following equations: 

1 ^ M h-Mk^ + ^ + V^U/4 , 

yr X a -Bm(ki + </>)-y/=AU/A _ yr X a — Xfi — 

which are coupled nonlinear equations called Lieb-Wu equations [21]. The Lieb-Wu equations 
have not been solved analytically. However it predicts some important results on thermody- 
namic properties of the system through Takahashi's string hypothesis [16, 28-30]. In terms of 
the solutions {h,X a } of the Lieb-Wu equations (2.3), energy eigenvalues are written as 

H^\k,X; S f NtM = E\k,X;s)% tM , E = -2 £ +</>) + \u{L-2N). (2.4) 

l<i<N 
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Recall that the Bethe states (2.2) give only the eigenstates with ^ 2M < N ^ L. In order 
to construct other eigenstates, we need to consider the symmetries of the system. 
Hereafter we shall sometimes abbreviate the superscript 4> m \k,\;s)^ N M . 

2.2 Symmetries 

The ^/-independent symmetries of the Hubbard model are classified in [7]. First we review 
the symmetries connected with the spin and charge degrees of freedom [26, 27, 32, 33]. Define 
the following operators related to the spin degrees of freedom: 

s z ■= j E ( n n- n n)> S + : = E C *W> S - : = ( s +) j - ( 2 - 5 ) 

They give a representation of the algebra sl(2) on the Fock space V. Since all the operators 
{S Z ,S±} commute with (2.1) for an arbitrary value of cf), it is said that the system has 
spin-s/(2) symmetry. One finds another representation of sl(2) related to the charge degrees 
of freedom, 

^■■=\ E (l-r^-nn), 7?+ := £ e^ 2 ^%c ih r?_ := ( V+ )l (2.6) 

Note that all the operators in (2.6) are commutative with {S z , S±} (2.5). It is easy to see that 
the operator rj z also commutes with for an arbitrary value of </>. However other operators 
rj± commute with only for the special values of 4> satisfying i + jSZ. Thus the system 
has charge-s£(2) symmetry if ^ + £ Z, and it reduces to charge-'u(l) symmetry given by 
rj z for other values of <p. Combining the above two kinds of sl{2) symmetries, we see that 
the system has so(4)(~ sl(2) © sl(2)) symmetry if ^ + — £ Z, and u(V) © sl(2) symmetry 
otherwise. For simplicity, we call the former so(4)-case and the latter sZ(2)-case. For the later 
discussion, we define the Casimir operators for each sl(2), 

V 2 ■■= \ (v+V- + V-V+) + Vl S 2 := \ + SS + ) + S 2 Z , 

which are employed to see the dimension of each representation. 
Next we introduce the following three operators: 

Ts-.= n p iV ,n, T Ph -.= n n(4+ c -)' n n p ^-v» 

where Pis-jt ■= 1 — (c\ s — cj- t )(cj s — Cjt) and [^J denotes the greatest integer in x. One notices 
that T s induces a spin-reversal transformation, T p h a particle-hole transformation and T r a 
reflection of the lattice. Direct calculation shows 

rp—1 TT rp TJ rp—1 TT rp TT rp — 1 TT rp TJ 
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For the system with even L, i.e., a bipartite lattice, we also introduce 

Combining these operations, we obtain the following transformation properties of the Hamil- 
tonian (2.1): 

T-^Ts = H+, T^T^H^T T T ph = H^, for even and odd L, 
T^T-^H^TrTpt = for even L. (2.7) 

One notices that for both even and odd L, the system has spin-reversal symmetry. While the 
system with even L has particle- hole symmetry brought by the transformation TbT r T p h, the 
system with odd L does not have particle-hole symmetry for generic values of <f> except the 
special values satisfying ^ + y 6Z 

2.3 Construction of non-Bethe states 

We construct eigenstates that are not included in the Bethe states (2.2). First we consider 
the system with so(4) symmetry, i.e., both spin-sZ(2) and charge-s/(2) symmetries. Here 
we recall that the AB-flux parameter <\> takes a special value: ^ + — G Z for even or odd 
L. Then, as shown in [5] for even L with <p = 0, the Bethe states (2.2) characterized by 
finite-valued solutions of the Lieb-Wu equations (2.3) are so(4)-highest, i.e., rj + \k, A; s)jv,m = 
S + \k, A; s)at,m = 0. Since 

r} 2 \k, A; s) n ,m = rj(rj + l)|fc, A; s)n,m, S 2 \k, A; s) n ,m = S(S + l)\k, A; s) n ,m, (2.8) 

with i] = ^(L—N) and S = ^(N—2M), the Bethe state \k, A; s)tv,m is the highest weight vector 
of an (L — N + 1)(N — 2M + l)-dimensional highest weight so(4)-module. Here we note that 
rj + S is an integer for even L [33], while it is a half-integer for odd L. The so(4)-descendant 
states of the Bethe state |fe, A; s)n,m, 

(i]-) n (S-) m \k,\;s) NM , (0<n^L-N,0<m^N-2M), (2.9) 

are also eigenstates of (2.1) (see Fig. 1). They are energy degenerate with |fe, A; s)n,m- It 
is easy to see that such application of the lowering operators ry_ and S- to the Bethe states 
produces other eigenstates than those with ^ 2M ^ N ^ L. Essler, Korepin and Schoutens 
counted the number of Bethe states (2.2) and their so(4)-descendant states (2.9) under the 
string hypothesis to show the combinatorial completeness of Bethe ansatz [3]. 

Next we consider the system only with the charge-«(l) and the spin-s/(2) symmetries. The 
Bethe states (2.2) are s/(2)-highest and satisfy only the second relation in (2.8), which means 
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— o(3.2) 
S- 

(5,2) • -0(5,3) 

(7,3) O -0(7 ; 4) 

so (4) -case sl(2) -case 

Figure 1: Bethe states and non-Bethe states for L = 5. The Bethe state with (N, M) = (3,1) 
denoted by a closed circle produces five so(4)-descendants denoted by open circles if the system 
has so(4) symmetry. As the so(4) symmetry reduces to sl(2) symmetry, the six eigenstates 
form three doublets of si (2). Then we need one more Bethe state with (5,2). Note that the 
eigenstate with (7, 3) is not a Bethe states and is constructed through the transformation 
T T T ph . 

that the Bethe state \k,\;s)fj M is the highest weight vector of an (N — 2M+l)-dimensional 
highest weight sZ(2)-module. The sZ(2)-descendant states of the Bethe states 

(S„) m \k,\;s)% iM , (0 <m^N -2M), (2.10) 

are eigenstates of (2.1) (see Fig. 1). One notices here that the Bethe states (2.2) and 
their sZ(2)-descendants (2.10) do not produce the eigenstates with L < N ^ 2L since the 
application of the lowering operator 5_ does not change the number of electrons. Such 
eigenstates with L < N ^ 2L are constructed as follows: i) by applying the transformation 
TbT r T p h to the Bethe state \k, A; s)2 L _ N £_m> (L < N ^ 2M ^ 2L) if L is even, or by applying 
T T T ph to the Bethe state \k, \; s)^l N>L _ M , ( L < N < 2M < 2L ) obtained by the Lieb-Wu 
equations (2.3) with 4> + tt instead of <f> if L is odd, we get the lowest weight vector of a highest 
weight sZ(2)-module; ii) the application of the raising operators (S+) n , (0 < n ^ 2M — N) to 
the lowest weight vector produces other degenerate eigenstates. 

Even if we are interested only in the system with = 0, the Bethe ansatz method needs 
the system with <fi = it for odd L. One of the main purposes in this article is therefore to 
discuss whether or not the above procedure produces all the eigenstates. 

3 Lieb-Wu equations with one down-spin 

By employing a graphical approach, we exactly count the number of finite-valued solutions 
for the Lieb-Wu equations (2.3) in the case when the system contains an arbitrary number of 
up-spins and one down-spin, i.e., N ^ 2 and M = 1. We remark that such exact analysis is 
presented in [2, 24] for the case = 0. In this section, we assume ^ <p < x since has the 
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sm(ki+(t>) -A = jCot(^), exp (v^4 ^ A;^) = 1. (3.2) 



same energy spectra as those of H^^n. In the case M = 1, the string hypothesis [28] predicts 
that two types of solutions exist; one is the solution with only real wavenumbers {fcj} and 
another includes two complex wavenumbers. We investigate such types of solutions below. 

3.1 Real k solutions 

First we consider the real solutions. For M = 1, the Lieb-Wu equations (2.3) reduce to 
e ^ = A-sin( W )-^/4 (i = 12 ^ 

A - s in(fc t + 0)- yp[U/A _ l 

These are equivalent to the following equations: 

U fhL 

) , exp \y/-l 

We investigate the real solutions for the first equation 

sin(g+0)-A = jcot(^). (3.3) 
In the interval ^ q < 2ir, its right hand side has L branches 

Note that, with a given £ satisfying (3.4), the first equations in (3.1) are rewritten as 

, r- « „ f X — sin(fcj + 6) \ 
k { L = 2tt£-2 arctan [ — -^J , 

which are convenient to relate the £ to the (half-)integers appearing in the string hypoth- 
esis [28]. By regarding A as a real parameter, we seek a solution q of (3.3) in one of the 
branches (3.4). From the graphical discussion (Figure 2), the solution in the branch I is 
uniquely determined for arbitrary A under the following condition: 

d . . , d (U fqLw UL 

-1= mm — sin g + <p) - X) > max — — cot — = - — — . 

0^<2tt dq 0^q<2n dq\4 V 2 // 8 

Hence, for U > f , the solution of (3.3) can be written as an increasing function of A, i.e., 
q = qt(X). Given a non-repeating set {£i\l ^ i ^ N} C {^-^|1 ^ j ^ L} of the branches, 
the second equation in (3.2) is satisfied when ^ J2i Qe^X) G Z. The behaviour of the solution 
q = qe(X) tells us that 

l<i<N l<i<N 
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Figure 2: The generic behaviour of sm(q+<j)) — A and ^-cot(^) in (3.3) under the condition 
UL > 8 for L = 6. Horizontal dashed lines divide six branches as specified by (3.4). Each 
branch has an intersection which leads to a solution of (3.3). 

Thus there exist N — 1 values of A giving the following integer values for ^- ^ i q^.{\): 



2 

l<i<N 



1,2,...,A-U. 



Note that such {A} and integers {m} are in one-to-one correspondence due to — > 0. 
As a consequence, the solutions {ki = q£.(\),\} are specified by the set of indices {£i,m}. 
The number of possible {£i,m} is given by (\r)(N — 1). Here we note that the number is 
consistent with the formula Z(L;N,M) in [3,28]. 

Proposition 3.1. For U > f , the Lieb-Wu equations (2.3) with M = 1 have (^)(N - 1) 
real solutions [2, 24] ■ 

3.2 A>A-string solutions 

Next we consider the solutions including a couple of complex wavenumbers. We assume the 
form of solutions as 

ki £ M/2^M, (i = 1, 2, . . . , N - 2), k N -i = C- V^Te, k N = C+ V^l£, 

where ^ £ < 27r and £ > 0. Note that /cat_i and /cat form a complex conjugate pair which 
is referred to as fe-A-2-string. Then the first set of equations in (3.1) are rewritten in terms 
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of real variables as follows 



U /k-L\ 

sm(ki+(f)) - A = —cot ( —5— ) , (i = 1,2,... ,iV- 2), (3.6a) 



4 V 2 

sin(C + ^)coshg- A= ^ w,T! (CL) 777v ( 3 - 6b ) 

4 cosh(£L) - cos(C-L) 

cos(C+^) sinhe = u fn (gL) ^ n - (3.6c) 

4 cosh(4L) — cos((X) 

On the other hand, the second equation in (3.1) is equivalent to the following condition: 

2ir 

ki + 2( = —m, with m = 0, 1, . . . , NL - 1. (3.7) 

In the same way as the previous case of section 3.2, if we consider a solution of each equation 
(3.6a) in one of the branches (3.4), and the solution can be written as a function of A. Given 
a set ^ i ^ N — 2} of non-repeating indices specifying the branches (3.4), we express 
the solutions of (3.6a) as k{ = q£.(X), (1 ^ i ^ N — 2). Then, from the relation (3.7), the £ is 
also written as a function of A, 

C = C(A):=^m-i £ qUX ), (3.8) 

for fixed {£i} and m. 

For an illustration, we consider (3.6b) and (3.6c) in the case N = 2. Since ( does not 
depend on A in the case, the equations (3.6b) and (3.6c) decouple into the following: 

A = sin(^m+0)cosh£, sinh£ = -- JZ__/(2)(£), (3.9) 

\L J 4cos(-^m+0) 

where 

(2) sinh(^L) _ fcoth(£L/2) for m £ 2Z, 

^ '~ cosh(^L) - (-l) m ~ jtanh(£L/2) for m € 2Z + 1. 

We seek a solution of the second equation in (3.9) through graphical discussion. Since 
f^ 2 \0 > when £ > 0, we need the condition ^ < ^-m + (f) < ^ so that the second 
equation of (3.9) has a solution. If such m is even, it is straightforward that the second 
equation in (3.9) determines a unique solution £(> 0) since lim^oo f( 2 \£) = 1. For odd m, 
the equation has a unique solution if the condition 

1 dj^Q . / £7 # (2) /n ^ UL 

1 = ^^ (0)< -SSi (-Mf^jlT (0) ) = T 

is satisfied. The number of allowed values of m here depends on L and <j), 

\{h -^ + i|i = i,2,...,L-i}, ^-^ez, 
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Figure 3: The behaviour of sinh£ and ~ 4cos (^ + ^) /a(0 in (3.10a) under the conditions UL > 8 
and < C + ^ < X" r ^ ne £' an< ^ £" res P ec tively denote the maximum and the turning point 

4cos«"+0) * 



° f ~ 4msf^l /°(0 111 tlle CaSeS < « < 1 



(2) (2) 

Let £m denote the solution specified by the m. By substituting the solution £m into the 
first equation in (3.9), one immediately obtains A. Recall that, in the cases of 4 — —(f) G Z, 
i.e., even L, (respectively, odd L) and <fi = or ^, (respectively, = w or |j), the system 
has the so(4) symmetry. Thus, as the so(4) symmetry reduces to the spin-s/(2) symmetry 
through the change of AB-flux strength <f>, the number of solutions for the Lieb-Wu equations 
increases. 

Let us consider the case N > 2. By inserting (3.8) into (3.6c), we have 

SinK = -4cos(aA)+0) /a(e) ' (3 - 10a) 

/.(£):= , « = (-rcos(^ £ %i (A)). (3.10b) 

cosh(£L)-a V2 !^_ 2 1 

Note that \a\ ^ 1. One sees that, for a fixed a, the function / a (£) has the following properties: 
i) /o(0) = 0, ^r(O) > lim^oo /„(£) = 1; ii) if a < 0, / a (£) is monotonically increasing and 
concave; iii) if o > 0, f a (0 has a single positive maximum at £'(> 0) and a single turning point 
at £"(> CO- These properties are sufficient to discuss the solution £ of (3.10) for arbitrary 
A. From the graphical discussion similar to the case N = 2 (see Figure 3), this determines a 
unique £ as a function of A under the condition 

1 = ;t — (0) < mm — ~. 77 — rr-rrlO) = mm tt; 7 = —^Ti 

d£ yj kKi V 4cos(C+</>) d£ v V |a|<;i4(l-a) 8 

if and only if § < C(^) + 4> = | £V <?4(A) < ^f-. By using (3.5), it is sufficient to have 
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a unique solution for (3.10) that the integer m satisfies 



l<j<Af-2 l<i<N-2 



that is, 



' { Ei + 5) + f - ^ + i|i = 1, 2, . . . , L-iv+i}, § - ^ g z, 
{ Ei (4 + 5) + li - i^J + j\j = 1, 2, • • • , L-iv+2}, § - fa £ z. 



(3.11) 



(2) 

Note that liniA^ioo £(A) = £ „ ,,, ,i v which is well-defined for the above m. We see that, 
for the values of m given in (3.11), the equation (3.6b) with £(A) and £(A) 

A = sin(fm+<M£,. % .(A))cosh fe(A)) - 77 g^A^i^) 

U 2 2- «. U 4 cos (| % (A)) - cosh (C(A)L) 

=:<?({%(A)U(A)), (3.12) 

determines A. In fact, since <?^(A) and £(A) are continuous functions of A and the function g 
satisfies the following: 

(2) 

m-Ei(^±§) 



lim </({%(A)U(A)) =g ({^(4±|)U 



Here 5 is a continuous and finite function with respect to A. Hence there exists a solution A 
for the equation (3.12). 

Proposition 3.2. For U > f , the Lieb-Wu equations (2.3) with M = 1 have (^^(L-N+l) 
k-A-2-string solutions if the system has the so(A) symmetry, and they have ( N L _ 2 ){L — N + 2) 
k -A- 2- solutions, otherwise. 

One notices that, only for the system with the so(4) symmetry, the number of k-A-2- 
solutions is consistent with the string hypothesis [28]. Note that, for < U < some of 
the A;-A-2-strings may disappear. In Appendix A, we numerically investigate the case N = 2 
and show that, for < U < the fc-A-2-strings with odd m disappear, while additional 
real solutions appear [5]. For the system with only the spin-sl(2) symmetry, the Lieb-Wu 
equations have more fc-A-2-solutions than those expected by the string hypothesis. 

3.3 Completeness of Bethe ansatz for L = 3 

Applying the above results, we now show that all the eigenstates can be constructed through 
the Bethe ansatz method in a simple case. We consider the case L = 3 and <p 7^ |^ when 
the system does not have so(4) symmetry but the spin-sl(2) symmetry. We note that the 
completeness of eigenstates in this situation has not been discussed in the literature. 
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The number of Bethe states \k, A; s) N M , (0 ^ 2M ^ N < 3) is exactly calculated as 
follows: the case M = is trivial since the eigenstates are those of lattice free fermion 
system; for the cases (N, M) = (2, 1) and (3, 1), we have obtained the following formulas from 
Proposition 3.1 and 3.2: 



jj (Bethe states) 



3 \ (N-l 
Nj V 1 

3 \ fb-N 
N-2J \ 1 



for real solutions, 



for fe-A-2-string solutions. 



Since each Bethe state | k, A; s)^ M corresponds to the highest weight vector of a highest weight 
s£(2)-module, we should count their s/(2)-descendant states 

(S_) n |M;< )Af , (0<n^N-2M). 

The eigenstates with 4 ^ N ^ 6, which are not Bethe states nor their s/(2)-descendant 
states, are constructed through the transformation T r T p h as we have described in Section 2.3. 
Indeed, by applying T r T p h to the Bethe state \k, A; s)q-n 3-Mi (4 ^ N ^ 2M ^ 6) of the 
system described by the Hamiltonian H^+k, we get the eigenstate T r Tph|fc, A; s)g^J^- 3 _ M of 
with 4 ^ N ^ 2M ^ 6 which is the lowest weight vector of a highest weight sZ(2)-module. 
We also count the eigenstates 

(S + ) n T r T ph \k, A; s)£^ 3 _ M , (° < n < 2M " iV )- 

Table 3.3 indeed shows that we obtain 64 = 4 3 = dimV eigenstates, which give a complete 
system of the Fock space V. 



4 Combinatorial completeness of Bethe ansatz 

The Bethe ansatz method was first introduced in the case of one-dimensional spin-^ isotropic 
Heisenberg spin chain [1]. Bethe assumed the string hypothesis and estimated the number 
Z(N; M) of solutions for the Bethe equations with M down-spins on an iV-site chain as 



v n (™) 



(4.1) 



{M n } m^l 



where M n denotes the number of n-strings composing a solution and P n = N — 2M + 
^Em(>n)( m ~~ n)M m . He obtained the following summation formula: 



Z(N;M) = ( 



N 
M 



N 
M - 1 



(4.2) 
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AT 

iv 


M 


P. AT 

— iv 


o — ikz 


type of solutions 


fH-t>etrie) 


si\Z ) sym. 


tl (state) 












real 


1 


1 


1 


1 









real 


3 


2 


6 


2 









real 


3 


3 


9 


2 


1 






real 


3 


1 


3 


2 


1 






fc-A-2-string 


3 


1 


3 


3 









real 


1 


4 


4 


3 


1 






real 


2 


2 


4 


3 


1 






rv-1 I-Z-Slllllg 


u 


9 


1 9 


4 


2 


2 


1 


real 


3 


1 


3 


4 


2 


2 


1 


fc-A-2-string 


3 


1 


3 


4 


3 


2 





real 


3 


3 


9 


5 


3 


1 





real 


3 


2 


6 


6 


3 








real 


1 


1 


1 












64 



Table 1: Enumeration of eigenstates for L = 3 and </> / |J . 

which implies that the number Z(N; M) of Bethe states is interpreted as the multiplicity of 
(N— 2M+l)-dimensional irreducible s/(2)-modules in the tensor product of N two-dimensional 
irreducible sZ(2)-modules. By taking into account that the Bethe states are s£(2)-highest and 
generate (N — 2M + 1) s/(2)-descendant states, the completeness of Bethe ansatz is shown in 
a combinatorial way as 

^2 (N -2M +l)Z(N;M) =2 N . 

0^M^[N/2] 

where [x\ denotes the greatest integer in x. 

It is known that, in general, solutions of Bethe equations do not have the nature assumed 
in the string hypothesis [4,8,9]. Indeed, for N > 21, some of the 2-string solutions are 
redistributed to real solutions that are not counted in Z(N; M) [4]. Hence, when we actually 
employ the the formula Z(N; M) to show the completeness of Bethe ansatz, we must regard 
such redistributed real solutions as 2-string solutions. 

We apply the techniques developed in [12, 13] to the Hubbard model with so(4) symme- 
try. Indeed a new proof for Essler-Korepin-Schoutens' combinatorial completeness of Bethe 
ansatz [3] is obtained as a corollary of the relation among the characters of so(4)-modules. 
Moreover, based on the results in the previous section, we propose the conjectural formula 
related to the combinatorial completeness of Bethe ansatz for the system with only the charge- 
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u(l) and spin-sZ(2) symmetry. In both cases, we obtain the formulas corresponding to (4.2), 
which has not been established even for the so(4)-case. 

4.1 Kirillov's power series and Q-system 

First we give three lemmas introduced in the case of one-dimensional isotropic Heisenberg 
spin chain [12, 13]. Detailed proofs are given in [12, 13]. Let a±, 0,2, ■ ■ ■ , tty be a set of integers 
for I » L. Define a set of formal power series {<f n (z n , z n+ i, . . . , zi)\l ^ n ^ 1} by 

Mz) := (l-z)- a - +2M ~\ 

ip n (z n ,z n+1 , ...,z t ):= t/j n (z n )ip n+1 ((l- z n )~ 2 z n+1 , . . . , (l-z n )~ 2 ^~ n hi). 
Lemma 4.1. The power series (p n (z n , . . . , z\) has the following expression: 

Mzn,..., Zl )= £ II ( Vmia ^ + Mm )z^...z^, for l<n<Z, 

where 

V n {a n ) = a n -2M+2 (m-n)M m . 

m>n 

Proof. The case n = I is given by the formula, 



z m . 



m^O 

Then the case 1 ^ n < I is proved by induction on n. □ 
Introduce the variables {z^\l ^n^l,0^k^n} through 

4 0) = *., 4* ) = (i-^" 1) )- 2(B -* ) 4*- 1) , ^ io^n. 

(k) 

Lemma 4.2. The power series <pi(zi, . . . ,z\) is rewritten in terms of the variables {zn }, 

Mzi,..., zi )= n (l-z^r^™- 1 - 

Proof. Since 

Mz { r\- ■ ■ , = ^(^^(a-rv 2 ^. • • • . (i-^-^-^-^r 1 ^ 

= (1 _^-l) ) -a fc+ 2M- W(z W ) 

for 1 ^ ^ / — 1, the lemma is proved. □ 
Define the polynomials {Q n = Q n {t)\n G Z^o} through the recursion relation, 

Qn+2 = Qn+1 — tQn, Qo = Ql = 1- (4.3) 
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Lemma 4.3. 

->2 



i) Q n = Qn+iQn-i + t n , for n > 1, (4.4a) 

n) Mt,t 2 ,--.,t l )= [J (Qn + lQn 2 Qn-lT an+2M -\ (4.4b) 

(l-v) n ~ 1 (l-v n+1 ) V 

Hi) Q n (t(v))= y - {1 _ v2) n Wkere : = (l +t ,)2 - ( 44C ) 



Proof, i) Use induction on n. 

ii) Set = z n = t n , (1 ^ n ^ /) in {zi^}. One obtains the following relations: 

i-zt 1] =Q k+ iQlQk-i, 4*> = Qll { r k) QT ~ k ~ 1] t n , for 10<n. 

Combining these with Lemma 4.2, we can prove (4.4b). 

hi) One can directly verify that the Q n (t(v)) satisfies the recursion relation (4.3) with t = 
t(v). □ 

The relations (4.4a) are called the Q-system of type sl(2), which is a key object in [12]. 
Indeed the expression (4.4b) produces an identity among the characters of sZ(2)-modules. 
The Q-system also plays a significant role in the combinatorial identities associated with the 
XXZ-Heisenberg spin chain and its generalizations [14, 17-19]. 

4.2 Combinatorial formulas 

Using the above lemmas, we discuss the Hubbard-case. In the similar way, we define (p' n (z n , z n +±, 
and V n (a' n ) by replacing a n and 2M with a' n and N in ip n (z n , z n+ i, ... ,zi) and V n (a n ), re- 
spectively. Define 

<p(s, t) := (1 + s) Vi( A s 4 t 2 , s 2 V)<^(t, t 2 ,..., t l ). 
The following is straightforward from Lemma 4.1: 

=Y J ( L )n ( v ' n ^ a ' n ^ M ' n \( Vn ^ +Mn \ Nr+2 ^^ 



{M n ,M{i) 



Then the coefficient of s N t M , (0 ^ 2M ^ N < L) in tp(s,t) is expressed by 



Z(L,{a n ,a' n };N,M):= £ (tf ) II 



V' n {a' n )+M' n \ (V n {a n )+M n 



M' n J\ M n 



N=N r +2Y. nM' n 
M=J2 n (M n + M' n ) 



where the sum runs over all configurations {iV r , M n , M' n ^ 0} such that N = N r +2 Y^ n >o n ^-n 
and M = ^ n>1 ^(M n + M' n ). We calculate explicit forms for the power series ip(s,t) after 
taking the special values of {a n } and {a' n }. 
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Introduce 



{iV r ,M n ,A<} v ' n^l 

iV=JV r +2E ™A#n 
M=E n(M„+M4) 



{JV r ,M n ,A<} v ' n^l 
JV=JV r +2E 
M=£ n(M„ + i\<) 



where 



P' n = L-N+2 {m-n)M' m , P n = N-2M+2 £ (m-n)M m . 

m>n m>n 

Note that P' n , P n ^ due to < 2M < JV < L. 

The N, M) (4.5a) is the very number of Bethe states for the Hubbard model estimated 
under the string hypothesis [3,28,30]. In terms of the string hypothesis, N r denotes the 
number of real fc's, M n the number of A-ra-strings, and M' n the number of fc-A-2n-strings. We 
have verified in the previous section that, if the system has so(4) symmetry, the number 

gives the correct number of solutions for the Lieb-Wu equations (2.3) with M = 1. In fact 
the first and second terms in the above Z(L; N, 1) correspond to the following two cases: N 
real fc's (N r = N) with a A-l-string (Mi = 1), and N - 2 real fc's (N r = N — 2) with a 
fc-A-2-string (M[ = 1), respectively. 

We now propose the Z(L; N, M) (4.5b) as a formula counting the number of Bethe states 
for the system with charge-u(l) and spin-s/(2) symmetries. Indeed 

is consistent with the number and the string-type of solutions for the Lieb-Wu equations (2.3) 
with M = 1 in the sZ(2)-case. The first term corresponds to the case of A" real fc's (N r = N) 
with a A-l-string (Mi = 1), and the second term to the case of Af — 2 real fc's (N r = N — 2) 
with a fc-A-2-string (M[ = 1). We note that, for L < N < 2L, we interpret Z(L;N,M) as 
the number of the lowest weight vectors T r T p h|fc, A; s)2~l-n l-m- ^° derive Z(L; N, M) (4.5b) 
for M ^ 2 from Takahashi's string center equations [28, 30], we need to appropriately extend 
the region of the allowed (half-)integers characterizing the Bethe states. 

In both the so(4)- and s/(2)-cases, one must also take a redistribution phenomenon into 
consideration [3, 4, 8, 9]; what it means here will be more clear in Appendix A. 
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Proposition 4.4. We have the following identities: 



)u N v M , 



i) (l+u) L {l+uv) L (l-u 2 v){l-v) = ^ Z{L;N,M)i 
i) {l+u) L {l+uv) L {l-v)= Z{L;N,M)u N v M . (4.6) 



Proof. By using Lemma 4.3, we have 

V(s,t) = (l + S) L [J (QUl^Qn-J^^-HQn+lQ^Qn-l)- "^- 1 , 



where := Q n (s t). Through the change of variables 

u(l + «) t; (l-tAQq-t;) 

1 + (l + t>) 2 ' (l + «V) 2 (l + i;) 2 ' 

the coefficient of s N t M in <p(s,t) is calculated as 
Z(L,{a„,a;};iV,M) = s= Res =o 

- Res / ( 1 + n ) L ( 1 + ^) L TT iqi n'~ 2 0' ,Y< +N - l (0 ^Q~ 2 1 r a " +2M_1 

M)T, | X2M+2 (l-n 2 t;)(l-^) \ rf" ^_ 

(i + t ,)iv+i 1 (i + ^) 2 (i + t>) 2 ;^+ l « M + l 

- Res / (1 + ^) L (1 + H L (1-^)(1-^) -TT fl ^.n^^ n )tn \ ^ ^_ 

- n= ^ S =0 l (l + « 2 ^ (1 + V) N j 1 J JuN+lvM+l' 

where we have introduced 

&n = ~a n + 2a n _i - a n _ 2 , b' n = -a' n + 2<4_! - a' n _ 2 , (1 < n ^ /), 

with a_i = ao = a'_ 1 = a' = 0. Note that, in the third equality, we have employed the 
assumption I >>> L. Setting a n j»i = L and a^>x = N, we have —6^ = b 2 = L, —b\ = b 2 = N 
and b' n>3 = b' n>3 = 0, which proves the identity i) in (4.6). And, setting a n j>i = L + n and 
°n>i = we have b[ = —L — 1, b' 2 = L, —b\ = b 2 = N and b' n ^ 3 = b' n ^ 3 = 0, which proves 
the identity ii). □ 

Through the change of variables u = xy^ 1 and v = y 2 , we have 

(x+x- 1 +y+y- 1 ) L {x-x- 1 )(y-y- 1 )= Y N, M)x~ L+N ~ 1 y~ N+2M ~ 1 , (4.7a) 

(x+x -i +y+ y-i ) i ( y_ y -i) = £ Z(L;iV,M)x- L+ V JV+2M " 1 . (4.7b) 

0^JV^2L 
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First we consider the relation (4.7a). If we express the left-hand side of (4.7a) as F(x,y), 
we find the property F(x _1 ,y) = F(x,y _1 ) = —F(x,y). Then the first relation (4.7a) is 
rewritten as 

„L-N+l _„-L+N-l N-2M+1 _„.-N+2M-l 

( X+X -i+y + y-y= I! Z(L;N,M) X V - — * , (4.8) 

O^M^ \_N/2\ 

where [x\ denotes the greatest integer in x. Let be an (n+l)-dimensional irreducible sl{2)- 
module associated with the charge-^ (2) symmetry and V n that associated with the spin-s£(2) 
symmetry. We introduce an (n+l)(m+l)-dimensional irreducible so(4)-module by the tensor 
product Vn^m = <8> V m . Through the representation (2.5) and (2.6) of so(4), the Fock space 
V of the L-site system is isomorphic to the tensor product (V^o© Vb,i)® L as an so(4)-module. 
Note that 

r? > ) = i(? + 1 ) |w) ' S » = ?(? + 1 )l^' for \v) G K,m C V. 

We now decompose the Fock space V into the direct sum of V n ^ m . From the relations (2.8), 
each Bethe state \k, A; s) n,m corresponds to the highest weight vector belonging to Vl-n,n-2M = 
V2rj,2S- The characters of the so(4)-module V n , m are calculated as 

n+l _ — n— 1 „.m+l _ ..— m— 1 

i -rr x x y y 

cn v n m — — : — : . 

x — x L y — y 1 

To be precise, x = e Al and y = e Al where both A[ and Ai are the fundamental weight of 
sl(2) and they are orthogonal to each other. One notices that, in terms of the characters, the 
identity (4.8) can be rewritten as 

(chF li0 + chV ,i) L = Z(L;N,M)chV L —N,N—2M- (4.9) 

0<MsJ lN/2\ 

Theorem 4.5 (Multiplicity formula). The multiplicity of the irreducible component Vl~n,n-2M 
in the tensor product (V^o © Vq,i)® L is given by Z(L; N, M), 

(V lfi © F ,if L = Z(L;N,M)V L -N,N-2M- (4-10) 

Next we turn to the relation (4.7b). If we express the left-hand side of (4.7b) as F(x,y), 
we find F(x,y~ 1 ) = —F(x,y). This gives 

_ JV-2M+1 -JV+2M-1 

(l+x( J /+y- 1 )+a: 2 ) L = £ Z(L; iV, M) x N . (4.11) 

0<M< \_N/2\ 
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Let c V be a subspace of the Fock space F with iV electrons and let vi N) C y (A ° be 

an (m + l)-dimensional irreducible s/(2)-module related to the spin-sZ(2) symmetry in V^ N \ 
The VL N) can be regarded as a 1 x (m+l)-dimensional irreducible (u(l) © sZ(2))-module by 
employing the representation (2.5) of sl(2) and taking u(l) = C(L — 2i] z ) with (2.6). Here 

(L-2 Vz )\v) = N\v), S 2 \v) = ^(^ + l)\v), for \v) G C V. 

The Fock space V of the L-site Hubbard model is isomorphic to the tensor product (Vq ^ © 
as a (u(l) © sl(2)) -module. We decompose the Fock space F into the direct 
sum of From the first relation in (2.8), each Bethe state \k, X;s)n,m is the highest 

weight vector belonging to V^_} 2M = V^jp ■ The characters of Vm^ are calculated as 

m+l _ -m-1 

47^=1^ ^ . 

In terms of the characters, the identity (4.11) can be rewritten as 

(chy (0) + chy 1 (1) + chy (2) ) i = £ z(L ; iv,M)ch^ ) 2M . (4.12) 

Theorem 4.6 (Multiplicity formula). T/ie multiplicity of the irreducible component V^} 2 m 
in the tensor product (V (0) © © V O (2) ) 0L zs ^roen fry Z(L; N, M), 

(y (0) © © vi 2) r L = © N, M)V^} 2M . (4.13) 

0^JV^2L 
O^M^ lN/2\ 

Corollary 4.7 (Combinatorial completeness). We have 

i) dimy= {L—N+1)(N—2M+1)Z(L; N, M), 

0<M< lN/2\ 

%%) dimy= ^ (N-2M+l)Z(L;N,M). (4.14) 

0^JV^2L 
0^M< \_N/2\ 

Proof. Consider the limit x,y — > 1 in the identities (4.8) and (4.11). □ 

The identity i) in Corollary 4.7 reproduces the combinatorial completeness of Bethe states 
for the Hubbard model with so(4) symmetry obtained by Essler, Korepin and Schoutens [3]. 
The factor (L— N +l)(N — 2M+1) in i) corresponds to the dimension of the highest weight 
so(4)-module Vl-n,n-2M with the highest weight vector \k, A; s)n,m- 

In Essler-Korepin-Schoutens' proof of i) in Corollary 4.7, they take the sum on N after 
taking that on M. In our proof, the sums on N and M are taken "simultaneously" in the 
level of characters. 
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The factor (N—2M+1) of the identity ii) in Corollary 4.7 for ^ N ^ L is the dimension of 
the highest weight (u(l) © s/(2))-module V^} 2M with the highest weight vector \k, A; s)n,m- 
For L < N ^ 2L, the factor (N— 2M+1) should be interpreted as the dimension of the highest 
weight (n(l) © sZ(2))-module V^} 2M with the lowest weight vector T r T p h|A;, A; s}2L-N l-m- 
If even L, the identity ii) can be rewritten as 

dimV = 2 1 ~ 6n ' l (N-2M+1)Z(L;N,M), 

0^M< \_N/2\ 

by considering the particle-hole symmetry of the system (2.7). Thus we speculate that the 
identity ii) in Corollary 4.7 accounts for the combinatorial completeness of Bethe states for 
the system with the charge-u(l) and the spin-s£(2) symmetries. 

The identities (4.7) also enabled us to get the explicit form of Z(L; N, M) through the 
binomial theorem. 

Corollary 4.8. We obtain the summation formulas for Z(L; N, M) and Z(L; N, M), 

ii) Z { L-,N,An = (i)(,\,)-(J,)(, i,,V (415) 



MJ\N-MJ \M-lJ \N-M+1 

5 Summary and concluding remarks 

In the framework of Bethe ansatz, we have studied the Hubbard model with the AB-flux that 
controls the symmetry of the system. In Section 3 we have shown the existence of solutions 
for Lieb-Wu equations with an arbitrary number of up-spins and one down-spin. We have 
found that the number of fc-A-2-solutions increases as the so(4) symmetry reduces to the spin- 
sl(2) symmetry (Proposition 3.2). The number of Bethe states is consistent with the string 
hypothesis only in the so(4)-case. In Section 4 we have investigated the combinatorial formulas 
giving the combinatorial completeness of Bethe states. We have shown that the number of 
Bethe states can be interpreted as the multiplicity of irreducible components in the tensor 
products of so(4)-modules (Theorem 4.5). Essler-Korepin-Schoutens' combinatorial formula 
is reproduced by the relation (4.9) among the characters of so(4)-modules (Corollary 4.7). An 
advantage of our approach is that we can obtain the summation formula (4.15) for Z(L; N, M). 
We have also proposed a new combinatorial formula derived from the relation (4.12) among the 
characters of (u(l)©s£(2))-modules. The formula is related to the combinatorial completeness 
of Bethe ansatz in the s/(2)-case (Corollary 4.7). It should be remarked that, in Section 3, 
we have not proved the uniqueness of solutions. Although the Lieb-Wu equations may have 
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other solutions that we have not expected, the combinatorial formulas introduced in Section 4 
give an evidence that, for M = 1, we have found solutions enough to verify the combinatorial 
completeness of Bethe ansatz. The problem is open for M ^ 2. 

The combinatorial completeness of Bethe ansatz has not been discussed for the one- 
dimensional isotropic Heisenberg spin chain with twisted boundary conditions. Kirillov's 
identity [12] also produces the following formula: 

e ii( p " + m" 

where P n = N — 2M + < 2Yl m (>n)(' m ~ n )M m - It is expected that, if the redistribution 
phenomenon of solutions for Bethe equations [4, 8, 9] is taken into consideration, the formula 
corresponds to the system with twisted boundary conditions. We remark that the formula 
also appears in the different context [20]. 
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A Redistribution phenomenon 

In Section 3, we exactly show the existence of solutions of Lieb-Wu equations with M = 1 for 
U > -|. There, real solutions have been specified by non-repeating indices {£i\i = 1,2,..., N} 
and m. But, for < U < -|, real solutions with repeating indices may appear at the same 
time as /c-A-2-string solutions disappear. Such redistribution of type of solutions is observed 
in the isotropic Heisenberg model for a large number of sites [4,8,9]. Here we numerically 
investigate such phenomenon for the Hubbard model with L = 20 and N = 2 [3] . 

As we have already mentioned, fc-A-2-string solutions with odd m may disappear for 
< U < -| (see Figure 3). For each m, the critical value of U is exactly given by [3], 




Plotted on Figure 4 are the center ( = j^m, (m = 11, 12, ... , 29) of fc-A-2-string solutions and 
redistributed real solutions {k\, k 2 } for L = 20 and N = 2 in varying the value of U. Plotted 
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on Figure 5 are their imaginary parts. The /c-A-2-string solutions with odd m disappear for 
U < [/M and, at the same time, real solutions with repeating indices {£1, £2; m} = {y, tt; m} 
appear. Note that, as U — > 0, all the fc-A-2-string solutions on Figure 4 approach to the 
wavenumbers of lattice free fermion system. Thus, when we apply the combinatorial formulas 
in Corollary 4.7 to the case, we must count the number of Bethe states by regarding the real 
solutions with repeating indices as /c-A-string solutions. 



2k 



8/L 



Figure 4: The real parts of solutions for Lieb-Wu equations with L = 20 and N = 2. The 
dashed lines correspond to the centers £ of fc-A-2-string solutions with even m, and the dots 
express the centers £ of fc-A-2-string solutions with odd m and their redistribute real solutions. 
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Figure 5: The imaginary parts £ of solutions for Lieb-Wu equations with L = 20 and N = 2. 
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